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Abstract 

Tne  applicability  and  domains  of  validity  for  the  Bom  and  the  Ryiov 
methods  in  scattering  theory  are  established,  with  mathematical  rigor,  by 
comparing  successive  terms  of  the  Bom  and  the  Rytov  series  calculated  for 
wave  propagation  in  a  homogeneous  dielectric  half-space  and  a  homogeneous 
dielectric  slab,  for  which  the  permittivities  are  assumed  to  have  a  low  contrast 
over  that  of  free-space.  While  the  (first-order)  Rytov  approximation  is  superior 
to  the  (first-order)  Bom  approximation  when  it  is  applied  to  estimate  the 
scattered  field  in  the  homogeneous  dielectric  half-space,  both  approximations 
are  inapplicable  to  estimate  the  scattered  fields  when  the  dielectric  slab  is  very 
thick. 
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1 .  Introduction 


There  is  much  interest  in  modelling  the  propagation  and  scattering  of 
waves  in  inhomogeneous  media,  to  solve  both  direct  and  inverse  scattering 
problems.  The  physical  properties  of  these  media  can  be  extracted  or  predicted 
from  the  measurable  quantities  such  as  the  scattered  fields  and  intensities  which 
can  be  derived  from  the  models.  Usually,  approximations  are  sought  in  order  to 
make  the  scattered  fields  (which  can  be  expressed  as  the  Fredholm  integral 
equations  of  the  first  kind)  more  tractable  for  numerical  computation.  Two 
first-order  approximations  (i.e.,  the  Bom  and  the  Rytov  approximations)  are 
most  commonly  exploited  in  the  Bom  and  the  Rytov  methods.  However,  the 
application  of  the  Rytov  method  [1-68]  has  raised  considerable  debate  over  its 
relative  merit  and  applicability  over  the  Bom  method  and  vice-versa  [46].  For  a 
scatterer  of  compact  support  (with  d  as  the  size  of  its  largest  dimension),  the 
empirical  remark  is  usually  made  that  the  Bom  approximation  is  valid  only 
when  the  scatterer  is  small  on  the  scale  of  the  incident  wavelength  X  in  free- 
space  and  the  strength  of  its  scattering  function  V  is  weak,  i.e.,  low  contrast 
with  respect  to  free-space.  This  can  be  summarized  by  requiring  that  Vkod  «  1 
where  ko  =  InlX.  On  the  other  hand,  the  criterion  for  the  Rytov  approximation 
to  be  applicable  requires  that  the  fluctuations  in  V  be  slow  on  the  scale  of  \  but 
the  strength  of  V  is  not  necessarily  small  or  low  contrast.  In  other  words,  this 
physical  interpretation  of  the  validity  of  the  Rytov  approximation  is  based  on 
the  requirement  that  the  absolute  value  of  the  rate  of  change  of  the  complex 
phase  of  the  field  is  ver\’  small  compared  with  koV^. 

In  order  to  refine  and  clarify  the  domains  of  validity  for  both  methods 

and  thus  establish,  more  precisely,  the  situations  for  which  method  is  preferable 

over  the  other,  we  consider  two  deterministic  scattering  cases  in  this  paper.  Tlie 
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first  scattering  geometry',  described  in  section  2,  consists  of  free  space  to  one 
side  and  a  dielectric  !  ilf-space  of  permittvity  on  the  other.  We  model 

the  propagation  of  the  field  across  the  interface  between  the  free-space  and  the 
dielectric  half-space  and,  then,  make  the  comparison  on  the  results  obtained 
from  the  Bom  and  the  Ryiov  methods.  In  section  3,  we  present  and  compare  the 
results  derived  from  both  methods  for  the  second  scattering  geometry  which 
consists  of  a  dielectric  slab  of  finite  extent  in  the  propagation  direction.  Based 
on  these  comparisons,  we  are  able  to  delineate  the  appropriate  domains  of 
validity  for  both  methods  in  section  4  and  provide  concrete  conclusions,  which 
can  be  generalized  to  other  scatterers,  in  particular,  scatterers  of  compact 
support. 

2.  One-dimensional  Homogeneous  Dielectric  Half-space 

Consider  a  time-harmonic  plane  wave  with  radial  frequency 

CO  incident  from  free-space  (region  0  for  z  <  0)  into  Hr-,  one-dimensional  (1-D) 
homogeneous  dielectric  half-space  (region  1  for  z  >  0)  [see  Figure  1].  The 
wavenumber  in  free-space  is  kg  (=  co'y/  p.o£o  )  where  |J.o  and  Eq  are  the 

permeability  and  the  permittivity  for  free-space,  respectively.  Assume  that  the 
dielectric  half-space  has  the  permeability  Pi  and  the  permittivity  eo(l  +  V) 

where  the  scattering  function  V  is  real  and  positive. 

2.1.  Exact  Solutions  derived  from  the  Boundary  Conditions 

The  close-form  solutions  for  the  total  fields  'Fo(z)  and  T'i(z)  in  free-space 
(region  0)  and  in  the  dielectric  half-space  (region  1)  can  be  easily  derived  from 
the  boundary  conditions: 


and 


%(0)  =  %{0) 


(1) 


1  d%(z) 


^0  dz 


z  =  0 


1  d4^,(z) 


dz 


z  =  0. 


(2) 


These  boundary  conditions  can  be  applied  for  the  acoustic  case  or  the  TE  case  in 
the  electromagnetic  wave  theory.  For  Pi  =  Po»  die  wavenumber  in  the  dielectric 

half-space  becomes  kj  =■  koVl+V  .  After  matching  the  boundary  conditions, 
we  obtain 


%(Z)  =  +  R(k,) 


and 


z<0  (3) 

z  >  0.  (4) 

where  the  reflection  and  the  transmission  coefficients  R(kj)  and  T(ki)  are  given 


W,(z)  =  T(k0 


below: 
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T^(^'>  =  k„.k 
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When  V  «  1,  we  can  expand  ^0(2^)  ^nd  T'jCz)  into  Taylor  series  in  terms  of  V, 
i.e., 
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z  >  0.  (8) 


The  above  two  power  series  are  exactly  equal  to  the  Bom  for  the  Neumann) 
series  which  can  also  be  derived  in  scattering  theory'  under  the  Bom  method, 
shown  belov. . 


2.2.  The  Bom  Method  for  the  1-D  Dielectric  Half-space  Case 

In  scattering  theory,  To(z)  and  'Fi(z)  satisfy  the  scalar  wave  equations 


,  +  kj  %(2)  -  0: 
dz 

z  <  0 

(9) 

and 

d^'Fifz)  , 

2  +ki'Fi(z)  =  0; 

dz 

z  >  0. 

(10) 

Since  ki  is  equal  to  koV  1  +  V  ,  Eq.  (10)  can  be  rewritten  as 

'j  +q'i',(z)--ko‘V'P,(z); 
dz 

z  >  0. 

(11) 

where  the  term  on  the  right-hand  side  of  Eq.  (11)  is  the  so-caUed  secondary 
source  due  to  the  permittivity  contrast-EoV  of  the  dielectric  half-space  relative 

to  that  of  free-space.  Using  the  1-D  unbounded  scalar  Green's  function  Gfz,  z') 

which  is  governed  by  the  equation 

d‘Gfz.  z') 

^2  +  ko  G(z,  z  )  -  5(z  -  z  ), 

dz 

(12) 

we  c'm.  express  ’FpCz)  and  4^i(z)  as  an  inhomogeneous  Fredholm  integral 
equations  cf  ihe  first  kind,  namely. 


and 


T„(z)  =  J  dz'G(z,  z-)VM',(z');  2<0  (13) 

oo 

4^,(z)  =  4^1%  -  ^  I  dz'Giz,  2')V4^i(z');  z  >  0  (14) 
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where  the  unperturbed  fields  4^0  \z)  and  \z)  are  both  equal  to  e  .  By 
substituting  'Fi(z)  of  Eq.  (4)  into  the  integrals  of  Eqs.  (13)  and  (14)  and  using 
the  1-D  unbounded  scalar  Green’s  function  [69], 


G(z,  z')  = 


ikolz-z'I 


2ko 


(15) 


it  can  be  shown  that  %(z)  and  'Fi(z)  derived  from  integral  representations  of 
Eqs.  (13)  and  (14)  are  consistent  with  that  of  Eqs.  (3)  and  (4)  derived  from 
boundary  conditions  of  Eqs.  (1)  and  (2).  Note  that,  in  Eq.  (14),  the  integral 
with  limits  (0  to  «>)  is  divided  into  two  integrals  with  limits  (0  to  z)  and  (z  to  '^) 
corresponding  to  two  regions  {z'  <  z)  and  (z'  >  z)  as  indicated  by  the  absolute 
value  of  Iz  -  z'\  in  Eq.  (15).  The  principle  of  limiting  absorption  is  used  for  all 
derivations  to  discard  the  oscillating  terms  of  the  form  when  z  — > 
[25,41].  The  Bom  series  of  Eqs.  (7)  and  (8)  in  terms  of  the  small  parameter  V 
can  be  derived  from  the  perturbation  method  [2.6]: 
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and 
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=  H-,  (2)  +  %,  (r)  +  'F„  (z)  +  4'„  (z)  +  ...; 


z>0  (17) 


nBA  nBA 

where  yqs  (z)  and  'Pu  (z)  (n  =  1,  2,  ...)  are  the  n-th  order  Bom  approximated 
scattered  fields  in  regions  0  and  1,  respectively.  With  straightforward  but 
tedious  calculation  up  to  the  third-order  Bom  approximation,  we  have 


=  (t  -  i  V^oZ  -  V^o 


z<0,  (18a) 


z<0,  (18b) 


z  <  0,  (18c) 


z>0,  (19a) 


z>0,  (19b) 


5V^  .  5V^ ,  3V^ , 2  ' 
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-  i  ^ko  z>0.  (19c) 


Eqs.  (18a)-(19c)  are  exactly  the  same  as  the  corresponding  terms  in  Eqs.  (7) 

oo  nBA 

and  (8).  Under  the  constraint  V  «  1,  the  Bom  series  X%s  (^z)  is,  but  not  the 

n=l 
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Bom  series  (z)  when  VkoZ  >  1,  uniformly  valid  for  all  z.  For  a  fixed  V 

n=l 

“  nBA 

and  ko,  this  nonuniformity  problem  occurred  in  the  Bom  series 

n=l 

imposes  a  severe  limitation  on  the  propagation  distance  z  for  the  Bom  method 
to  be  useful. 


2.3.  TTie  Rytov  Method  for  the  One-dimensional  Dielectric  Half-space  Case 


Let  MoCz)  and  %(z)  be  defined  as 


(0)  i(Pn(z) 

%(z)  =  %  \z)  e  •, 


z<0  (20) 


z>0  (21) 


where  (Pq(z)  and  (Pj(z)  are  the  complex  phase  functions  for  regions  0  and  1, 

respectively.  For  the  case  that  kj  =  koV  1  +  V  7  exact  solutions  of  Eqs.  (3)  and 
(4)  yield  exact  expressions  for  9q(z)  and  9j(z),  namely, 


9o(z)  =  -  i 


,  r'Fo(z)'|  r  ,^-i2koz' 

Ln  -“Tot —  =  - 1  Ld  1  +  R(k,)e  ; 

1%  U)J 


z  <  0  (22) 


%iz) 


cp,(z)  =  -  i  Ln  =  -  i  Ln[T(k,)]  (k,  -  ko)z;  z  >  0  (23) 

[%  (z)J 

where  Ln[X]  is  the  natural  logarithm  of  X.  If  V  «  1,  we  can  expand  (Po(z)  and 
(Pj(z)  into  Taylor  series  in  terms  of  V,  i.e., 

,  ,  .  rv  5vh  -ilkoz  .i4koZ  .  V-'  -i6koZ 
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z>0.  (25) 


Under  the  definitions  of  Eqs.  (20)  and  (21),  we  rearrange  the  wave  equations  in 
Eqs.  (9)  and  (10)  and  obtain  [35,64] 

d‘r  (0)  1  7  (0)  (0) 

(z)(p„(z)J  +  T„  (z)<p„(z)  =  -  %  (z);  z  <  0  (26) 


d‘"r.„(0)  1  ?  (0)  7  (0)  rdcp,(z)]  (0) 

(z)(p,(z)J  +  T,  (z)<p,(z)  =  i  k^  VH-,  (z)  -  i  %  (z); 

z  >  0,  (27) 

Therefore,  with  the  application  of  the  1-D  unbounded  scalar  Green’s  function 
G(z,  z'),  we  can  express  tp^Cz)  and  (Pj(z)  as 
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It  is  easy  to  show  that  the  exact  solutions  of  <Pq(z)  and  9,(z)  in  Eqs.  (22)  and 

(23)  can  also  be  derived  from  Eqs.  (28)  and  (29)  by  using  Eq.  (15)  and  the 
exact  phase  changes; 


and 


d(Pp(z)  2koR(ki) 

^  1  +  R(k,) 

d(p,(z) 
dz  ^ 


z  <0  (30) 


z>0.  (31) 


Eqs.  (24)  and  (25)  can  also  be  derived  from  the  Rytov-approximated  solutions 
of  Eqs.  (28)  and  (29),  i.e., 


^0  (z) 


dz'G(z,z') 


d9'^'^z')  dtp^j^'^z')  dtp'J^'tz') 
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=  (Pj  (z)  +  tpj  (z)  +  (pj  (z)  +  ...; 


z>0  (33) 


nRA  nRA 

where  (Pq  (z)  and  (p  j  (z)  are  the  n-th  order  Rylov  approximated  complex 

phase  functions  in  regions  0  and  1,  respectively.  After  lengthy  but  simple 
calculations,  we  obtain 
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The  factor  of  two  in  Eqs.  (34c)  and  (35c)  is  essential  for  the  derivation.  With 

increasinsly  cumbersome  formulations,  one  can  follow  the  same  iteration 
scheme  to  derive  higher-order  terms  of  the  Rytov  series  for  tp^lz)  and  (p^(z). 

Since  the  penurbation  expansions  for  (pjz)  in  Eqs.  (32)  and  (33)  are  based  on 
the  criteria  that  (z)l  «  (z)l  (m  =  0,  1  and  n  =  1,2,  ...)  which,  in 

turn,  leads  to  the  constraint  V  «  1,  it  imposes  no  restriction  on  the  propagation 
distance  z  for  the  Rytov  method. 

From  Eqs.  (34a)-(35c),  we  have  shown  that  the  conventional  way  of 
omitting  [d(Po(z)/dz]^  and  [d(Pi(z)/dz]"  in  Eqs.  (24)  and  (25)  under  the  Rytov 

approximation  [1,2,6]  is  equivalent  to  neglect  the  multiple  scattering  processes 

[11]  which  correspond  to  higher-order  terms  (n  >  1)  in  Eqs.  (32)  and  (33). 
Also,  we  have  shown  that  the  omission  of  [d(pQ(z)/dz]^  and  [d(Pj(z)/dz]  in  the 

Rytov  approximation  has  nothing  to  do  with  the  magnitudes  of  the  complex 
phase  functions  as  long  as  the  contributions  from  [d(pQ(z)/dz]  and  [d(Pj(z)/dz] 

are  negligible  compared  with  koV  [Ref.  17,  Appendix  A].  TTierefore,  the 
argument  [7.29,32]  that  the  Rytov  approximation  is  not  valid  when  IcPq(z)I  ~  1 

and  l(Pj(z)l  ~  1  is  incorrect. 


2.4.  Relationship  between  the  Bom  and  the  Rytov  Methods 


W’e  have  recovered,  from  Eqs.  (18a)-(19c)  and  (34a)-(35c),  the  w'ell- 
known  relations  between  the  Bom-approximated  scattered  fields  and  the  Rytov- 
approximated  complex  phase  functions  [25-27],  namely. 


/  X  ^  IRA  .  . 


(36a) 


(36b) 


:eA  (0)  .  :ra  1 

^ms  (Z)  =  T^  (Z)  - - +1(P^  (Z)[, 


3BA  (0)  IR 

T„3  (z)  =  (z)!*--  "3,  -  cp'; 


^z)9-J^V)  +  i(p^j^^z)(,  (36c) 


and  so  on  for  higher-order  terms  in  region  m  (0  and  1).  In  order  to  compare 
domains  of  validity  for  both  Bom  and  Rytov  methods,  we  examine  the 
relationship  between  the  Bom  and  the  Rytov  solutions  shown  below  in  Eqs. 
(37a)-(38c)  where,  up  to  the  third-order  power  of  V,  the  scattered  fields 
derived  by  the  Rytov  method  can  be  reduced  to  that  of  derived  by  the  Bom 


method,  i.e.. 


(Z)^%  (z)Le 


.  IRA.  , 

(z) 


'{'i®*(z)  +  ofv'e 


2e*i3koz^ 


z  <  0,  (37a) 


z  <  0,  (37b) 


fl+2+3)RA 


T  IRA,  -  :RA,  ,  3RA,  .n 
x[%  (z)  +  9  0  ^  0 


(z)  =  %  (z)  ie 


z  <  0,  (37c) 


IBA  2BA  3BA  ,  -i3k/,zA 

=  (z)  +  4"os  (z)  +  (z)  +  G  ) 


IRA  ,  fO) 

4^1,  (z)  =  4>,  (z)Le 


.  IRA,  .  1 

l(p,  (z) 

-  1. 


(-  T+  i  +  o(v')  +  o(v'k^  z=) 

=  '{',‘f*(z)  +  o(v’)  +  o(v"k^  2^); 


z  >  0,  (38a) 


='(-  7+  >  TM  +  T-  >  Vm  -  TkJ  z’-)e‘'^  +  0(v')  +  o(v’kJ  z’-) 
=  'p)f*(z)  +  'fu*(z)  +  0(V^)  +  o(v\Jz2);  z>0,  (38b) 


f  r  IRA,  s  2RA,  ^  3RA,  si  ) 

=  il 

YlV  V  .  V  ,  v2  .  v2  2  5V^  .  SV^  3V3  2  2 

V-T+iTkoZ  +  — -i'4'koZ--g-koZ^--p-+i— koZ  +  — 

-  i  ^  kj  Z^)  +  OCV*)  +  o(v\^  Z‘) 

=  'Fi?^(z)  +  'Fu  ^(z)  +  ^(z)  +  oCV*)  +  o(yX  zO; 

In  the  Rytov  method,  the  scattered  fields  corresponding  to  the  Rytov  series 

n  iRA  n+2‘*-...+n)RA  (0)  f  r-  jRA,  ,1  -i  \  / 

Icpj^  (z)  IS  denoted  by  (z)  =  4"^  (z)  {  exp[il(p^  (z)]  -  1  )  (m  = 

J=1  Fl 

0,  1).  From  Eqs.  (37a)-(38c),  it  is  shown  that  the  multiple  scattering  up  to  the 

(l+2-^...-bn')RA 

n-th  order  is  included  in  the  Taylor  expansions  of  4^nis  '^hich  can  be 


n  iBA 

represented  by  the  sum  of  I'ifns  (z)  plus  an  infinite  number  of  redundant 

j=i 

higher-order  terms.  As  implied  by  Eqs.  (37a)  and  (38a),  it  is  incorrect  to 
address  the  Taylor  expansion  of  exp[pi(r)]  [Eq.  (12)  of  Ref.  13]  as  "a  multiple¬ 
scattering  interpretation  of  the  Rytov  solution"  because,  for  example,  the 
quantity  p]’(r)  =  [2k.o  jdv'G(r,r')N(r')Ao(r')/Ao(r)]",  which  is  the  square  of 

V 

pi(r)  representing  the  single  scattering  process,  does  not  correspond  to  the 

double  scattering  process  whereas  the  quantity 
4k^|dv'G(r,r')A^(r')  Jdv"G(r',r")A(r")Ao(r")Mo(r)  does.  Also,  from  Eqs. 

V  V 

(37a)-(37c),  we  have  shown  that  for  finite  n,  4^05^^^  '^"^^%)  is  inadequate  to 
approximate  the  scattered  field  'Ros(z)  (o*"  reflected  field  in  free-space) 
because  it  does  not  provide  the  correct  propagation  factor  t  ®  w'hereas  the 

n  jBA 

Bom  series  l4os  (z)  does  [as  stated  in  Ref.  41  for  n  =  1].  Hence,  althoush  both 

j=i 

methods  have  the  same  domain  of  validity  (defined  by  the  constraint  V  «  1) 
when  each  is  applied  to  estimate  'Fos(z),  the  Bom  method  is  superior  to  the 
Rytov  method.  On  the  other  hand,  to  estimate  the  scattered  field  'Ri5(z)  in  the 
dielectric  half-space,  the  Rytov  method  is  superior  to  the  Bom  method  because 

<»  oBA 

the  Bom  series  (z)  suffers  the  nonuniformity  in  convergence  if  VkoZ  >  1 , 

n=I 

w’hich  is  not  the  case  with  the  Rytov  method. 

(l+2+...+n)RA 

Note  that  the  redundant  terms  in  the  Taylor  expansions  of  %,s  (2:) 

(m  =  0,  1)  cancel  out  one  another  exactly  if  n  — >  Thus,  we  have 


3.  One-dimensional  Homogeneous  Dielectric  Slab 

In  this  case,  the  same  plane  wave  as  that  described  in  section  2  is  incident 
from  free-space  (region  0,  z  <  0)  into  the  homogeneous  dielectric  slab  (region 
1,  0  <  z  <  d)  with  permeability  po  permittivity  £(,(1  -f  V)  [see  Figure  2], 

The  boundary  conditions  at  z  =  0  and  z  =  d  are 

%{0)  =  4^i(0),  (40a) 


and 


4^,(d)  =  %(d),  (40b) 

6%(z)  d%(z) 

1  =  j 

02  2  =  0  dZ  2  =  0, 

d'F,(z)  d4^2(z) 

=-4“^ 

dz  z-d  dz  z  =  d 


where  ^2(2)  is  the  total  field  in  region  2  (free-space  for  z  >  d).  By  matching  the 
boundary  conditions,  the  total  fields  in  three  regions  are  given  as 


%(z)  =  e“^  + 


R(k,)(: 


-ikoz. 


D(kj) 


z  <  0, 


T(k,)^ikiz  R(ki)T(k,)e  ^  ^.ik,z 


D(k,) 


_  T(k,)[l-R(k,)]e‘*‘^'-^°’'‘ 

■  D(ki)  ^  ’ 


0<z<d,  (43) 


z>d,  (44) 


where  D(k])  h  1  -  [R(kj)]  .  To  reduce  the  total  fields  in  Eqs.  (42)  and 

(43)  for  the  dielectric-slab  case  to  that  in  Eqs.  (3)  and  (4)  for  the  dielectric  half¬ 
space  case  in  which  does  not  exist,  one  can  invoke  the  principle  of  limning 
absorption  [41]  to  assume  that  kj  contains  a  small  and  positive  imaginaiy  part 


such  that  '  =  0.  If  V  «  1  and  Vkod  «  1,  Eqs.  (42)-(44)  can  be  expanded 

into  the  Taylor  series: 

vPo(z)  {[- (,  .  e-od) 

.  i  (^  -  kod  ^  . ... )  z  <  o,  (45) 


Vl^«l  r  V  v2  5v3  fV^  5Vh^i2kocl  .v^  z  i2 

'^'i(z)  ~  LI -T+T-'sr+l^T^-  “6rje  ^  +1T6 Md  + 2^e 

v2  5Vh  f  v2  3Vh,2  2  .v\3  ikoZ 

rv  3v2  9v3  i2kod  .fW^  V^,  z, 

+  La  •  16  ■*■  64  +  64  ^  +  1  (^4  •  XjMd  -  2) 

,  2,j  Z.2l  _i2kod  ^-ikoz 

-XW*2)  J  e  °  e  °  +  ...;  0<z< 


i2kod 


0  <  z  <  d,  (46) 


V«1 

Vkod«l 


'^'z(z)  [1 .  (-  (1  -  .  i  (^ .  .  9m 

+  i  ^kod  +  (-  V  +  'W] -  i  ^ kj  d^  +  ...]  z  >  d.  (47) 


One  can  also  derive  Eqs.  (42)-(44)  by  using  the  inhomogeneous  Fredholm 
equations  of  the  first  kind  for  the  total  fields  in  region  0,  1,  and  2,  namely, 


u 

'i'Jz)  =  ’(z)  -  K  \  dz'GCz.  z^VH-.d'); 


(m  =  0,1,2)  (48) 


where  G(z,  z')  and  H^i(z)  are  given  in  Eqs.  (15)  and  (43),  respectively,  and 

(0)  (0)  (Ol 

(z)  =  (z)  =  Tn  (z)  =  e  .  With  tedious  but  straightforward  calculations, 

it  can  be  verified  that  successive  terms  in  the  Bom  series,  derived  iteratively  bv 

(0) 

replacing  %(z)  with  (z)  in  Eq.  (48),  are  exactly  equal  to  the  corresponding 
terms  in  the  Taylor-series  expansions  for  Eqs.  (42)-(44)  in  terms  of  V. 

Since  only  the  constraint  V  «  1  is  used  in  the  Bom  method  for  estimating 

op 

^m(z)  (m  =  0,  1,  2  ),  all  of  three  Bom  series  (z)  diverge  if  Vkod  >  1. 

n=l 

“o  nBA 

The  Bom  series  iT'u  (z)  is  subject  to  two  special  situations: 

n=I 

nBA  <»>  nBA 

(1)  if  z  =  0,  X'Eis  (z)  behaves  like  vT'os  (z)  and 

n=l  n=l 

^  nBA  "c  nBA 

(2)  if  z  =  d,  X^is  (z)  behaves  like  14^2$  (z)- 

n=l  n=l 

Therefore,  to  obtain  good  estimates  for  by  using  the  Bom  approximation 
requires  that  V  «  1  and  Vkod  «  1.  These  two  constraints  impose  a  severe 
limitation  on  the  domain  of  validity  for  the  Bom  method. 

Alternatively,  'Fo(z),  'Ei(z),  and  4^2(2)  can  be  defined  as: 

4'„(z)  =  (49) 

(0) 

where  4^^  (z)  and  <Pn,(^)  (m  =  0,  1.  2)  are  the  unperturbed  fields  and  the 

complex  phase  functions  for  regions  0  (z  <  0),  1  (0  <  z  <  d),  and  2  (z  >  d), 
respectively.  By  the  exact  solutions  of  Eqs  {42)-(44),  we  have 


9o(z)  =  -  i  Lrt  1  + 


R(ki)Vl  -  e 

D(ki) 


z<0.  (50) 


r  i2k,d 

|T(ki)  R(k,)T(k,)e  - 

(p,(z)  =  -  1  D(ki)  *  D(ki)  ^ 


+  fk,  -  ko)z;  0  <  z  <  d,  (51) 


rT(k,)r I  -  R(k,)]] 

(p,(z)  =  -  1  Lnj  +  (k,  -  ko)d; 


z>d.  (52) 


WTien  V  «  1  and  Vkod  «  1,  Eqs.  (50)-(52)  can  be  expanded  into  the  Taylor 
series  in  terms  of  V,  namely, 


V«l 

Vkod«l 


-  T^M  e““^] (i  -  e 


-i4koZ 


$ 


z<0,  (53) 


v«i 

Vkod«l 


Vkod«l  ^  3y2  /n  i2kod 

(pi(z)  ~  3r+  16+^  16+  le^odj^ 
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.  i  +  o(v')  +  o(v\?dO  +  o[\^(d  -  z?]  + ...; 


0  <  z  <  d,  (54) 


Vl^«I  yZ  y3  5y'‘^  .A'-  29  V 
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z  >  d.  (55 ) 


.  V\ ^i2kod 


.  \  idknd 
-ijl^e  +...: 


As  usual,  Eqs.  (50)-(52)  can  also  be  derived  from  the  following  integral 
representations,  i.e., 

0 

i  f  (0) 

(pjz)  =  -  -(oj—  dz'G(z,  z)  %  (z') 

%.(z)  J  L  dz  J 


if  ^  rd(p,(z')f  (0) 

+  Jdz'GCz,  z')|ko^V  -  (z') 

(z)q^ 

oo 

i  f  rd(p,(z')f  (0) 

-  “^(0) —  dz'G(z,  i')  ■  r.  %  feO;  (m  =  0,  1,  2)  (56) 

L  dz  J 


By  the  same  token  as  that  for  Eqs.  (34a)-(35c),  Eqs.  (53)-(55)  can  be  derived 
from  Eq,  (56)  by  iteration.  Hence,  as  in  Eqs.  (36a)-(36c),  the  same  well-knoum 
relations  between  the  Bom  and  the  Rytov  series  can  also  be  established  in  all 

oo  nRA 

three  regions.  From  Eqs.  (53)-(55),  it  is  shown  that  the  Rytov  series  X9i  (z) 

n=l 

is  also  subject  to  two  special  situations:  ‘ 

(1)  if  z  =  0,  LcD^'°^(z)  behaves  like  T(p"^'^(z)  and 

n=l'  n=l  ^ 

(2)  if  z  =  d,  icp"^^(z)  behaves  like  . 

n=l  rv=l  ‘ 

If  the  Rytov  approximation  is  able  to  provide  fair  estimates  for  (p^(z).  both 
constraints  that  V  «  1  and  Vkod  «  1  should  be  invoked.  Hence,  the  domain  of 


2  1 


validiiy  for  the  Ryiov  method  is  the  same  as  that  for  the  Bom  method  in  this 
dielectric-slab  case. 

It  is  worth  pointing  out  that  if  we  let  the  thickness  d  of  the  dielectric  slab 
approach  to  infinity  but  keep  V  and  kp  fixed,  we  are  unable  to  establish  the  same 

statement  as  that  in  the  dielectric  half-space  case  on  issues  of  applicability  and 
domains  of  validity  for  the  Bom  and  the  Rytov  methods.  This  is  attributed  to 

the  fact  that  we  have  invoked  the  principle  of  limiting  absorption  to  neglect  the 

highly  oscillating  terms  of  the  form  lim  which  should  exist  after  the 

2— 

integrations  are  carried  out  for  the  dielectric-slab  case,  even  when  d  is  very 
large  (but  finite). 

4.  Conclusion  and  Remarks 

We  have  shown  that,  to  approximate  the  scattered  field  (or  the  reflected 
field)  in  free-space  of  the  dielectric  half-space  case,  the  Bom  method  is  more 
suitable  than  the  Rytov  method  although  both  have  the  same  domain  of  validity. 
But,  to  compute  the  scattered  field  in  the  dielectric  half-space,  the  Rytov  method 

nBA 

IS  supenor  to  the  Bom  method  because  the  Bom  series  v'Rjj  (z)  suffers  the 

n=l 

nonuniformity  in  convergence  when  VkoZ  >  1.  In  the  dielectric-slab  case,  both 

Bom  and  Rytov  methods  have  the  same  (and  very'  narrow)  domain  of  validity. 
However,  when  Vk^d  >  1,  these  two  methods  are  inapplicable  because  the  Bom 

and  the  Rytov  series  diverge.  Despite  of  the  fact  that  the  principle  of  limiting 
absorption  provides  a  self-consistent  w’ay  of  deriving  the  exact  total  fields  for 
both  1-D  deterministic  cases,  it  is  unable  to  eliminate  the  secular  terms 
occurring  in  the  Rytov  series  for  the  dielectric-slab  case  when  d  becomes  very 
large  such  that  \kod  >  1.  In  summary,  the  domains  of  validity  for  the  Bom  and 


the  Rytov  approximations  are  described  in  tables  1  and  2  for  the  dielectric  half¬ 
space  and  the  dielectric-slab  cases,  respectively. 

The  above  statements  concerning  the  Bom  and  the  Rytov  methods  are 
also  valid  for  the  situation  when  the  dielectric  half-space  or  the  dielectric  slab 
has  a  permittivity  of  the  form  ei(l  +  V)  where  Ei  >  Eq-  In  this  case,  to  derive 

the  Bom  and  the  Rytov  series  from  inhomogeneous  Fredholm  equations  of  the 
first  kind  requires  a  set  of  appropriate  Green’s  functions  which  satisfy  the 
boundary  conditions  at  the  planar  interfaces.  We  have  shown  that  with  assistance 
of  these  sets  of  Green's  functions  for  1-D  deterministic  cases,  we  are  able  to 
develop  the  correct  integral  representations,  which  are  beneficial  to  both  direct 
and  inverse  scattering  theory,  for  the  exact  total  fields  in  different  regions  [70]. 
Note  that,  in  conjunction  with  the  application  of  the  statistical  theor}',  the  above 
theoretical  treatments  for  1  -D  deterministic  cases  can  be  extended  to  that  for  1  - 
D  stochastic  cases.  In  addition,  understanding  the  physical  interpretations 
inferred  from  analytical  expressions  for  these  1-D  cases  will  shed  some  light  on 
the  scattering  theory  for  3-D  cases  where  no  closed-form  expressions  exist. 
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Figure  Captions 


Figure  1.  Scattering  geometr\’  for  1-D  dielectric  half-space  case. 

Figure  2.  Scattering  geometry  for  1-D  dielectric-slab  case. 
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Table  1 


Domains  of  validity  of  the  Bom  and  the  Rytov  approximations 
in  the  dielectric  half-space  case  for  finite  ko  and  V  (V  «  1 ) 


Range 

App  rox  i  mat  i 

z  <  0 
for  all  z 

z  >  0 

z  «  (Vko)  ' 

z  >  0 

z  ~  (Vko)  ’ 

z  >  0 

z>(Vko)‘ 

Bom 

Y 

Y 

A 

N 

Rytov 

N 

Y 

Y 

Y 

A:  Ambiguous  to  estimate  the  scattered  field  or  the  complex  phase  function 
N:  Invalid  to  estimate  the  scattered  field  or  the  complex  phase  function 
Y;  Valid  to  estimate  the  scattered  field  or  the  complex  phase  function 


Table  2 


Domains  of  validity  of  the  Bom  and  the  Rytov  approximations 
in  the  dielectric-slab  case  for  finite  ko  and  V  (V  «  1 ) 


Range 

Approximation"*^".^ 

for  all  z 
d  «  (VV' 

for  all  z 
d  -  (Vk„)'‘ 

for  all  z 
d  >  W' 

Bom 

Y‘ 

A 

N 

Rytov 

Y 

A 

N 

3  2 


